Control of Hamiltonian chaos as a possible tool 
to control anomalous transport in fusion plasmas 
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It is shown that a relevant control of Hamiltonian chaos is possible through suitable small pertur- 
bations whose form can be explicitly computed. In particular, it is possible to control (reduce) the 
chaotic diffusion in the phase space of a Hamiltonian system with 1.5 degrees of freedom which mod- 
els the diffusion of charged test particles in a turbulent electric field across the confining magnetic 
field in controlled thermonuclear fusion devices. Though still far from practical applications, this 
result suggests that some strategy to control turbulent transport in magnetized plasmas, in partic- 
ular tokamaks, is conceivable. The robustness of the control is investigated in terms of a departure 
from the optimum magnitude, of a varying cut-off at large wave vectors, and of random errors on 
the phases of the modes. In all three cases, there is a significant region of maximum efficiency in 
the vicinity of the optimum control term. 
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I. INTRODUCTION 



Transport induced by chaotic motion is now a stan- 
dard framework to analyze the properties of numerous 
systems. Since chaos can be harmful in several contexts, 
during the last decade or so, much attention has been 
paid to the so-called topic of chaos control. Here the 
meaning of control is that one aims at reducing or sup- 
pressing chaos inducing a relevant change in the trans- 
port properties, by means of a small perturbation (either 
open-loop or closed-loop control of dissipative systems 
[lll3) so that the original structure of the system under 
investigation is substantially kept unaltered. Control of 
chaotic transport properties still remains an open issue 
with considerable applications. 

In the case of dissipative systems, an efficient strategy 
of control works by stabilizing unstable periodic orbits. 
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where the dynamics is eventually attracted. Similarly, 
a first idea to control Hamiltonian systems is to modify 
the parameters of the system in order to act on periodic 
orbits: One can enhance the stability of elliptic periodic 
orbits by zeroing their residues or by stabilizing hy- 
perbolic periodic orbits 0. 

Another idea to stabilize the system is to enlarge the 
phase space by coupling the system with an external sys- 
tem (and hence with additional degrees of freedom which 
makes the large system more reg ular) jl. These embed- 
ding techniques are similar to the above methods on the 
stabilization of unstable periodic orbits; they are based 
on the construction of a dissipative system from the orig- 
inal Hamiltonian system. The techniques developed for 
dissipative systems can thus be applied to this modified 
system, like for instance the targeting of periodic orbits. 
A different approach is to modify the Hamiltonian (or 
just the potential) to control the original system. This 
approach is useful when one is able to act on this sys- 
tem with an external forcing. The interesting point is 
that the Hamiltonian structure with its number of de- 
grees of freedom is preserved. So far, the modifications 
of the Hamiltonian that have been proposed in the liter- 
ature are: the modification of the integrable part of the 
Hamiltonian the control of a system with large and 
non-smooth external pulses , a localised control with a 
modification in some specific regions of phase space B , 
or a control using variations of the external field 0, ^3 ■ 
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However, we notice that most of the modifications of the 
potential that have been proposed so far are tailored to 
specific examples (with the exception of the optimal con- 
trol and the required modifications are large com- 
pared with the potential. 

Hamiltonian description of the microscopic origin of par- 
ticle transport usually involves a large number of parti- 
cles. Methods based on targeting and locking to islands 
of regular motions in a "chaotic sea" are of no practical 
use in control when dealing simultaneously with a large 
number of unknown trajectories. Therefore, the most ef- 
ficient procedure appears to be to control the transport 
process with a small perturbation, if any, making the sys- 
tem integrable or closer to integrable. In what follows we 
show that it is actually possible to control Hamiltonian 
chaos in this way by preserving the Hamiltonian struc- 
ture. We describe a general method for controlling nearly 
integrable Hamiltonian systems, and we apply this tec- 
nique to a model relevant to magnetized plasmas. 

Chaotic transport of particles advected by a turbulent 
electric field with a strong magnetic field is associated 
with Hamiltonian dynamical systems under the approx- 
imation of the guiding center motion due to E x B drift 
velocity. For an appropriate choice of the turbulent elec- 
tric field, it has been shown that the resulting diffusive 
transport is then found to agree with the experimental 
counterpart It is clear that such an analysis is only 
a first step in the investigation and understanding of tur- 
bulent plasma transport. The control of transport in 
magnetically confined plasmas is of major importance in 
the long way to achieve controlled thermonuclear fusion. 
Two major mechanisms have been proposed for such a 
turbulent transport: transport governed by the fiuctu- 
ations of the magnetic field and transport governed by 
fiuctuations of the electric field. There is presently a 
general consensus to consider, at low plasma pressure, 
that the latter mechanism agrees with experimental ev- 
idence _12J. In the area of transport of trace impurities, 
i.e. that are sufficiently diluted so as not to modify the 
electric field pattern, the E x B drift motion of test par- 
ticle should be the exact transport model. Even for this 
very restricted case, control of chaotic transport would 
be very relevant for the thermonuclear fusion program. 
The possibility of reducing and even suppressing chaos 
combined with the empirically found states of improved 
confinement in tokamaks, suggest to investigate the pos- 
sibility to devise a strategy of control of chaotic transport 
through some smart perturbations acting at the micro- 
scopic level of charged particle motions. 

As in the current literature the electric turbulent trans- 
port in plasmas is mainly addressed in the Eulerian 
(fluid) framework, let us first recall the difference be- 
tween Lagrangian and Eulerian descriptions of transport. 
We consider the advection of a scalar quantity 9{k, t) 
describing, e.g., the concentration of a passively trans- 
ported entity. In a given Eulerian velocity field v(x, t) 



the transport of 0(x, t) is described by 

^^1^ + v(x, t) ■ V^(x, t) = DV'9{^, t), (I.l) 

where D is a molecular diffusion coefficient. This equa- 
tion holds for both neutral fluids and plasmas. If the 
fleld v(x, t) is given indipendently from the field 0(x, i) 
Eq. (|I.1|I is linear in v(x, t). The complexity of the field 
0{x,t) will then depend on both the complexity of the 
field v(x, t) and on the molecular diffusion coefficient D. 
When considering the simulation of Eq. (|I.1|I . the mag- 
nitude of D will govern the mesh size to store the field 
6{x,t). For a vanishingly small diffusion D, there will 
be no cut-off of the small scales generated by the simu- 
lation. This will require an infinite storage capability to 
describe the complexity of ^(x, t) that can appear even 
for a relatively smooth velocity field. To evaluate this 
property, the most straightforward description is given 
by a Lagrangian approach. For the same transport pro- 
cess, the latter requires to solve the following equations 
of motion of a passive tracer (e.g. particle, ffuid drop) 
whose Eulerian concentration function is 6{x,t), 

x = v(x,t), (1.2) 

which in the case of a two dimensional incompressible 
Euler fiow can be given the form 

where tp denotes the stream function of the eulerian field 
v(x, t), the trajectory of the tracer is denoted by x(t) and 
V"*- = {—dy,dx)- What is remarkable here is the Hamil- 
tonian structure of the equations of motion ljl.3|) , where 
the stream function plays the role of the Hamiltonian 
function and x and y play the role of the canonically 
conjugate variables. With the exception of trivial veloc- 
ity fields V (like a uniform, parallel fiow) these equations 
of motion are in general nonlinear in the coordinates; 
in fact, if we even think of a simple vortex, we realize 
that Vx and Vy must contain at least one trigonometric 
function. Now, also without a standard (quadratic) ki- 
netic energy term, this kind of Hamiltonian dynamical 
system displays all the rich and complex phenomenology 
of the transition between regular and chaotic motions 
and between weak and strong chaos Thus, even 

in presence of rather regular Eulerian velocity patterns, 
the solutions of Eqs. (|I.2|) and HI.3|I can be very compli- 
cated, with apparently no relation left with v(x, f). In 
other words, chaotic Lagrangian diffusion can take place 
also in presence of rather simple Eulerian velocity pat- 
terns. For realistic simulation of Eq. (|I.1|I a finite mesh 
size must be introduced and accordingly the diffusion co- 
efficient D must reach a finite value to smear out the 
small scales that cannot be captured by the grid. If the 
velocity field is characterized by a large regular structure 
superimposed to small scale structures, the output of the 
simulation can lead to the absence of any diffusion but 
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the molecular one [3. The difficulty in the simulation 
of Eq. will then lead to an apparent conflict with 
a broad experimental evidence .15]. The most efficient 
means to address the transport of passive scalars in a 
given velocity field v(x, t) appears to follow a Lagrangian 
approach that allows one to describe the motion at all 
scales in space and time. The cost of this method will 
appear in the statistics that must be performed to ob- 
tain a general property of the system whenever a single 
trajectory does not allow one to capture the properties 
of all possible trajectories. 

When addressing plasma transport, Eulerian and La- 
grangian approaches are combined to provide an analy- 
sis of the transport properties. An equation similar to 
Eq. is coupled to a vorticity equation defining the 

field v(x, t) |J_6j. The Eulerian description is used to gen- 
erate the velocity field and a Lagrangian approach is used 
to follow trace impurities and trace Tritium that allows 
one to compare the simulations to experimental data |l7j . 

A close analogy exists between the equations of motion 
of passive tracers (|i.3ll and those of the guiding centers 
of charged particles moving in strongly magnetized plas- 
mas and in presence either of an electric field transverse 
to the magnetic field, or of an inhomogeneous component 
of the magnetic field itself. The electrostatic case J^] is 
modeled by 

■_ d fx\ _ c c f-dyV{x, y, t)\ 

where V is the electric potential, E = —WV, and B 
Bbz- The magnetic case is modeled by 

d f x\ vii 

^ Vjl(-dy%oi{x,y,t)\ 

RB\d^%oi{x,y,t) )' ^ - ' 

where V|| is the velocity along the field line, R the major 
radius of the torus and ^poi the poloidal magnetic flux 
divided by 2tt. In both cases, the physically remarkable 
phenomenon - in complete analogy with the Lagrangian 
diffusion of passive scalars - is that even in presence of 
rather regular space-time patterns of the electric fields 
or of the magnetic inhomogeneities, the charged parti- 
cles can diffuse across the magnetic field which ceases to 
be confining. The dynamical instability with respect to 
small variations of the initial conditions, known as deter- 
ministic chaos, is the very source of the enhanced cross- 
field diffusion; it is "intrinsically" non-coUisional and it 
turns out to be orders of magnitude larger than the coUi- 
sional one , sometimes even many orders of magnitude 
larger 

In this article, the problem we address is how to control 
chaotic diffusion in such Hamiltonian models. In some 
range of parameters, the problems can be considered as 
nearly integrable. We consider the class of Hamiltonian 
systems which can be written in the form H — Hq + eV 
that is an integrable Hamiltonian iJp (with action-angle 



variables) plus a small perturbation eV . 

The problem of control in Hamiltonian systems is the 
following one: For the perturbed Hamiltonian Hq + eV , 
the aim is to devise a control term / such that the dynam- 
ics of the controlled Hamiltonian Hq + eV + f has more 
regular trajectories (e.g. on invariant tori) or less diffu- 
sion than the uncontrolled one. Obviously / = —eV is 
a solution since the resulting Hamiltonian is integrable. 
However, it is a useless solution since the control is of 
the same magnitude of the perturbation. For practical 
purposes, the desired control term should be small (with 
respect to the perturbation eV) , localized in phase space 
(meaning that the subset of phase space where / is non- 
zero is finite or small enough), or / should be of a spe- 
cific shape (e.g. a sum of given Fourier modes, or with 
a certain regularity). Moreover, the control should be 
as simple as possible in view of future implementions in 
experiments. 

In Sec.ni we explain the control theory of nearly inte- 
grable Hamiltonian systems following Ref . ■ We show 
that it is possible to construct and compute analytically a 
control term / of order such that the controlled Hamil- 
tonian He = Hq + eV + f is integrable. In Sec. IIIII after 
defining the model of interest to our study in Sec. IIII Al 
we compute analytically the first terms of the expan- 
sion of the control term in Sec. IIII Bl Some properties 
of the control term are given in Sec. IIII CI A numerical 
study of the effect of the control term on the dynamics 
is done extensively in Sec. IIVI It is shown that the cho- 
sen control term is able to drastically reduce the chaotic 
transport. In Sec. EI '^e study the effect of some trun- 
cations that aim at either simplifying the control term 
or reducing the energy input to control the system: In 
particular, we show that reducing the control term to its 
main Fourier components or reducing the magnitude of 
the control term is sufficient to govern a significant de- 
crease of the chaotic transport. Thought, of course the 
optimal control is obtained with the full control term. 
These results indicates that this control of Hamiltonian 
systems is robust. 



II. CONTROL THEORY OF HAMILTONIAN 
SYSTEMS. 

In this section, following the framework of Ref. [23| we 
explain the control theory of Hamiltonian systems. Let 
A be the vector space of C°° real functions defined on the 
phase space. For H Q A, let {H} be the linear operator 
acting on A such that 

{H}H' = {H,H'}, 

for any H' E A, where {■ , ■} is the Poisson bracket. 
Hence ^ is a Lie algebra. The time-evolution of a func- 
tion V E A following the flow of H is given by 
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which is formally solved as 

V{t) = e*^^>F(0), 
if H is time independent, and where 



as 



Any element V ^ A such that {H^V = 0, is constant 
under the flow of H, i.e. 

Let us now consider a given Hamiltonian Hq G y^. The 
operator {Hq} is not invertible since a derivation has 
always a non-trivial kernel. For instance {Hq}Hq°' — 
for any a such that Hq"' e A. The vector space Ker {Hq} 
is the set of constants of motion. Hence we consider a 
pseudo-inverse of {Hq}- We define a linear operator T on 
^ such that 

{HoV r = {iJo}, (II.5) 

i.e. 

WVeA, {HoAHo.TV}} = {Ho,V}. 

The operator F is not unique. Any other choice F' satis- 
fies that the range Rg(F' — F) is included into the kernel 
Ker({Ho}2). 

We define the non-resonant operator J\f and the reso- 
nant operator TZ as 

= {Ho}T, 

where the operator 1 is the identity in the algebra of 
linear operators acting on A. We notice that Ea. (jII.5|) 
becomes 

{Ho}n - 0, 

which means that the range KgTZ of the operator TZ is 
included in KerjiJo}- A consequence is that TZV is con- 
stant under the flow of Hq, i.e. e*^^o>7^F = TZV. We 
notice that when {Hq} and F commute, TZ and Af are 
projectors, i.e. TZ^ = TZ and A/"^ = A/". Moreover, in this 
case we have RgTZ = KerjTJo}, i.e. the constant of mo- 
tion are the elements TZV where V € A. 

Let us now assume that Hq is integrable with action- 
angle variables (A, y>) G S x T' where B is an open set of 
and is the I dimensional torus. Thus Hq — Hq{A) 
and the Poisson bracket {H, H'} between two elements 
H and H' of A is 

SH H'\ = — ^ 
^ ' ^ dA dip dip OA' 

The operator {Hq} acts on V expanded as follows 



V 



J2 V^{A)e^-'^, 

kGZ' 



{Hq}V{A, 'p) = J2 • k Vi,iA)e'^-'^, 



where 



a;(A) 

A possible choice of F is 
TViA,ip)^ J2 



dHQ 
OA ■ 



T4(A) 
iijj(A) ■ k 



(n.6) 



We notice that this choice of F commutes with {Hq}. 

For a given V € A, TZV is the resonant part of V and 
J\fV is the non- resonant part: 

TZV^Y. ^k(A)x(u;(A) • k = 0)e'*"^, (IL7) 

k 

MV = J2 14(A)x(u;(A) • k ^ 0)e'''"^, (IL8) 



where x(q: = 0) vanishes when a =^ and it is equal to 1 
when a = 0. 

From these operators defined for the integrable part 
Hq, we construct a control term for the perturbed Hamil- 
tonian Hq + V where G .4, i.e. / is constructed such 
that Hq + V + f is canonically conjugate to Hq -\- TZV. 

Proposition 1 - For V A and F constructed from 
Hq, we have the following equation 



e^r^> {Hq + V + f) = Hq+ TZV, {11.9) 



where 



f{V) 



- ^-{rv} 



TZV 



l-e-{rn 
{TV} 



AfV-V. {11.10) 



We notice that the operator (1 — e ^'"^^)/{Fy} is well 
defined by the expansion 

Proof. Since e'^'"^^ is invertible, Eq. 1H.9(I gives 

f{V) = (e-{i^^> - 1)Hq + e-{i^^> W - V 

We notice that the operator e^'^'"^^ — 1 can be divided 
by {F^} 

fiV) = ^^^^^{TV}Hq + e-{i^^> W - V 

By using the relations 

{TV}Hq = {TV,Hq} = -{Hq}TV, 
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and 



we have 



{Ho}T=M, 



1 - p-{rv} 

/(V^)^e-{^^>W+ ^^^^ AfV~V. a 
The control term can be expanded in power series as 

n— 1 ^ 

We notice that if V is of order e, /(T^) is of order e^. 

Proposition 1 teUs that the addition of a well chosen 
control term / makes the Hamiltonian canonically con- 
jugate to Hq + TZV. It is also possible to show from 
Proposition 1 that the flow of Ho + V + f is conjugate to 
the flow of Ho + nV (see Ref. 



Proposition 2 



J{Ho + V+f} _ ^~{rV} t{Ho} MnV} JTV} 



The remarkable fact is that the flow of TZV commutes 
with the one of Hq, since {Ho}!^ — 0. This allows the 
splitting of the flow of H^ + TZV into a product. 

The notion of non-resonant Hamiltonian is defined by 
the following statement: 

Definition - Hq is non-resonant if and only if VA G 
B, u{A) • k = implies k = 0. 

If Hq is non-resonant then with the addition of a con- 
trol term /, the Hamiltonian Hq -|- y -|- / is canonically 
conjugate to the integrable Hamiltonian Hq -\- TZV since 
TZV is only a function of the actions [see Eq. HII.7|) ]. 

If Hq is resonant and TZV — 0, the controlled Hamil- 
tonian H = Hq + + / is conjugate to Hq. 
In the case TZV 



0, the series (IILII|I which gives the 
expansion of the control term /, can be written as 



or any function h{ijj ■ A), it is shown that if the frequency 
vector satisfies a Diophantine condition and if the per- 
turbation is sufficiently small and smooth, such a control 
term exists. An algorithm to compute it by recursion is 
provided by the proof. We notice that the resulting con- 
trol term N is of the same order as the perturbation, and 
has the following expansion 



^TZ{rV}V + Oie^), 



N{A) = TZV 



where we have seen from Eq. (III.7II that TZV is only a 
function of the actions in the non-resonant case. The as- 
sumption that <jj is non-resonant is a crucial hypothesis 
in Gallavotti's renormalization approach. Otherwise, a 
counter-term which only depends on the actions A can- 
not be found. 

Our approach makes possible the construction of a con- 
trol term in the resonant case. The controlled Hamilto- 
nian is conjugate to Hq -\- TZV where TZV depends on the 
angle and action variables in the resonant case. Therefore 
the controlled Hamiltonian is not integrable in general. 
The new term TZV which is always a conserved quan- 
tity is functionally independent of Hq since it depends 
on the angles. There exists a linear canonical transfor- 
mation (A', if') = (*TA, T^^ip) where T is a / x / matrix 
with integer coefficients and determinant 1 such that uj is 
mapped onto a new frequency vector which has its r last 
components equal to zero, where r denotes the dimension 
of {k G s.t. a; ■ k = 0}. In these new coordinates TZV 
depends only on r angles. This form of Hq-\-TZV is called 
the resonant normal form. The non resonant case occurs 
when r = 0. When r = I the normal form of Hq -f TZV 
depends only on one angle, so it is integrable. 
In what follows, we will apply the control theory to a 
resonant Hamiltonian which models the E x B drift in 
magnetized plasmas. 



III. CONTROL OF CHAOS IN A MODEL FOR 
E X B DRIFT IN MAGNETIZED PLASMAS 



/(F) 



(11.12) 



s=2 



where fs is of order and given by the recursion formula 

1 



-{TVJs-i}, 



(H.I3) 



where fi = V. 

Remark : A different approach of control has been de- 
veloped by G. Gallavotti in Ref. @. The idea is to find 
a control term (named counter term) depending only on 
the actions, i.e. to find TV such that 

H{A, if) - Hq{A) + F(A, v) - N{A), 

is integrable. For isochronous systems, that is 

i7o(A) =a;- A, 



A. The model 

In the guiding centre approximation, the equations 
of motion of a charged particle in presence of a strong 
toroidal mag netic field and of a time dependent electric 
field are n[ 



d 
di 



x=— =^E(x,<)xB = — 



<^/-dyVix,y,ty 

B\ d,V{x,y,t) 

(HI. 14) 

where V is the electric potential, E — —WV, and B = 
Be^. The spatial coordinates x and y where {x,y) G 
play the role of the canonically conjugate variables and 
the electric potential V{x, y, t) is the Hamiltonian of the 
problem. To define a model we choose 



kGZ2 



;(k)i 



, (HI. 15) 
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FIG. 1: Contour plot of V{x,y,t) given by Eg. am. 1611 for 

t = 0, a = 1 and TV = 25. 



where (/^k are random phases (uniformly distributed) and 
Vk decrease as a given function of |k|, in agreement with 
experimental data |2lj |. In principle one should use for 
w(k) the dispersion relation for electrostatic drift waves 
(which are thought to be responsible for the observed tur- 
bulence) with a frequency broadening for each k in order 
to model the experimentally observed spectrum 5'(k, w). 
In order to use a simplified model we use in this article 
ti;(k) = wq constant as a dispersion relation. The phases 
(/3k are chosen at random in order to mimic a turbulent 
field with the reasonable hope that the properties of the 
realization thus obtained are not significantly different 
from their average. In addition we take for |Vk| a power 
law in |k| to reproduce the spatial spectral characteristics 
of the experimental ^(k), see Ref. |21|. Thus we consider 
the following explicit form of the electric potential 



N 



1 



2tt J-^^^ [v? + m2)3/2 



n2 4-TTj,2< JV2 



27r 



{nx + my) + .^„,„ - ujot 



(III. 16) 

By rescaling space and time, we can always assume that 
L — 1 and = Stt. In what follows, we choose N — 25. 
Figure Q shows a visualization of the potential for t — {) 
and a = 1. Since the model is fluctuating in time, the 
eddies of Fig. ^ are rapidly modified in time and where 
a vortex was initially present, an open line appears, and 
so on. 

Two particular properties of the model, anisotropy and 
propagation have been observed : each image of the po- 
tential field shows an elongated structure of the eddies 
and superposing images obtained at different times a 
slight propagation in the y = x direction is found. How- 
ever, this propagation can easily be proved not to disturb 
the diffusive motion of the guiding centers. The prop- 
erty of propagation can be easily understood analytically. 
In fact, restricting ourselves to the most simplified case 
of an electric potential given only by a dominant mode 



(n = TO = 1) it is immediately evident that at any given 
time the maxima and minima of the sine are located on 
the lines y = ~x + constant. As the amplitudes are de- 
creasing functions of n and to, this structure is essentially 
preserved also in the case of many waves. The property 
of anisotropy is an effect of the random phases in pro- 
ducing eddies that are irregular in space. 

We notice that there are two typical time scales in the 
equations of motion: the drift characteristic time r^, in- 
versely proportional to the parameter a, and the period 
of oscillation of all the waves that enter the potential. 
The competition between these two time scales deter- 
mines what kind of diffusive behaviour is observed |llj| . 
In what follows we consider the case of weak or intermedi- 
ate chaotic dynamics (coexistence of ordered and chaotic 
trajectories) which corresponds to the quasi-linear diffu- 
sion regime (see Sec. lIV A^ . Whereas in the case of fully 
developed chaos that corresponds to the so-called Bohm 
diffusion regime one has to introduce a slightly more com- 
plicated approach (see remark at the end of Sec. IIII 



B. Computation of the control term 

We extend the phase space (x, y) into [x, y, E, r) where 
the new dynamical variable t evolves as T{t) = t + t(0) 
and E is its canonical conjugate. The autonomous 
Hamiltonian of the model is 



H{x,y,E,T)=E + V{x,y,T). 
The equations of motion are 



(III.17) 



dH 
dy 



1, 



dV_ . dH dV 

oy ox ox 

(III.18) 

and E is given by taking H constant along the trajec- 
tories. We absorb the constant c/B of Eq. (|III.14p in 
the amplitude a of Eq. ljIII.16|) . so that we can assume 
that a is small when B is large. Thus, for small values 
of a, Hamiltonian (|III.17|I is in the form H = Hq + eV, 
that is an integrable Hamiltonian Hq (with action-angle 
variables) plus a small perturbation eV. In our case 
Hq = E, i.e. independent of x,y,T, so that A = {E,x) 
and (f — (r, y) are action-angle coordinates for Hq {y can 
be considered as an angle but it is frozen by the flow of 
Hq). We could have exchanged the role of x and y. We 
have 



a;(A) = 



( OHq OHq 
\dE' dx 



(1,0), 



that is Hq is resonant [i.e. w(A) • k = does not imply 
k = ]. In order to construct the operators F, TZ and N 
we consider 

d__f d_ d_ 
dip \dT^ dy 



so 



{Hq} = UJ{A) 



d 



d 



dip dr 
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FIG. 2: Contour plot of /2 given by Eq. 1111.211 for a = 1 and 
N = 25. 



If we consider an element W{x,y,T) of the algebra A, 
periodic in time with period 1, we can write 



fcez 



and the action of F, TZ and TV operators on W is given 
by 



2ik'!T ' 



(III.19) 



and using the expressions of V and FV^, we have 

,2 1 



/2(X,' 



a 



E 



2(27r)3 „f^^ (n2 +m2)3/2(„2 ^^2)3/2 



{cos(27r(riix + miy) + (pmrm - ^ttt), 

sin(27r(n22; + TO22/) + 'Pn2m2 - Sttt)}, 

where {•, •} is the Poisson bracket for x, y coordinates, i.e. 
for two generic functions / and g depending on x,y,T, 
we have 

{/ 5} = — — - — — 

' dx dy dy dx 



From Eqs. (|III.16p and (|III.20p we obtain 

nini2 — Ti2mi 



f2ix,y,T) 8^ E^ (n2+m2)3/2(n2+TO2)3/2 

T12,TT12 

X sin[27r[(ni - ^2)2; + (rni - TO2)2/] + (ySmrm - 'j^'nams]- 

(III.21) 



We notice that for the particular model (|III.16p and for 
the particular choice of operator F given by Eq. (|III.19p , 
the partial control term /2 is independent of time. Fig- 
ure [2 depicts a contour plot of it. 
The computation of /a is given by 



nW^Woix,y), 
^fW = Y,Wk{x,y)e^'^''\ 



If we apply the operator F to y given by Eq. (|III.16|1 . 
we obtain 



FV" 



E 



1 



(27r)2 (^2+^2)3/2 

,i2 + „2<„2 

COS [2TT{nx + my) + (fnm — 2ttt] . (III. 20) 

Since V is periodic in time with zero mean value, we have 
TZV = 0. In this case, as we have seen in the previous 
section, Eqs. (jII.12|l - (jII.13|l give the expansion of the con- 
trol term. If we add the exact expression of the control 
term to Hq + V, the effect on the flow is the confinement 
of the motion, i.e. the fluctuations of the trajectories of 
the particles, around their initial positions, are uniformly 
bounded for any time |2Q|. 

In the present article we show that truncations of the 
exact control term / are able to regularize the dynamics 
and to slow down the diffusion. We compute the first 
terms of the series of the exact control term /2 and /s. 
From Eq. l|TrT3)l we have 

f2{V) ^ -^{TV^V}, 



and substituting the expressions ljIII.20p and (|III.21|I for 

TV and /2, one obtains 



h{x,y,T) 

E 

■n.\,mi ,712, m2 
Ti3,TrL3 



247r 

(niTO2 - min2) [(ni - ^2)^3 - (mi - rn2)"-3] 

{nl + m2)3/2(^2 ^ 77^2)3/2^^2 ^ fn'^)3/2 



X sin[27r(ni — 712 + n^)x + 27r(TOi — TO2 + ^3)1/ + 

fmrni - 'Vn^ra^ + ^n^ra^ " 27rT)] . (III. 22) 



The computation of the other terms of the series III. 121 



can be done recursively by using Eq. ljII.13|l (see also |2 
and the Appendix). 



C. Properties of the control term 

In this section, we first state that for a sufficiently 
small, the exact control term exists and is regular. The 
proofs of these propositions are given in the Appendix. 
Then we give estimates of the partial control terms in 
order to compare the relative sizes of the different terms 
with respect to the perturbation. 

Concerning the existence of the control term, we have 
the following proposition 
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FIG. 3: Contour plot of /a given by Eq. Jin.22t for t = 0, 

a = 1 and TV = 25. 



Another measure of the relative sizes of the control 
terms is by the electric energy density associated with 
each electric field V, /2 and /a. From the potential we 
get the electric field and hence the motion of the particles. 
We define an average energy density £ as 



1 

8^* 



E 



where E(a;,?/,f) = ~W. In terms of the particles, it 
corresponds to the mean value of the kinetic energy (i^ + 



y^) (up to a multiplicative constant), 
by Eq. (Himil . 



N 



£ = 



8^ ^ 



For V{x, y, t) given 



(III.23) 



Proposition 3 - If the amplitude a of the potential is 
sufficiently small, there exists a control term f given by 
the series fll.lS^) such that E + V + f is canonically con- 
jugate to E , where V is given by Eq. I^III.ldp . 

The proof is given in the Appendix. For = 25, it is 
shown that the control term exists for a < 7 x 10^^. As 
usual, such estimates are very conservative with respect 
to realistic values of a. In the numerical study, we con- 
sider values of a of order 1 . 

Concerning the regularity of the control term, we no- 
tice that each term fs in the series (|II.12|I is a trigono- 
metric polynomial with an increasing degree with s. The 
resulting control term is not smooth but its Fourier co- 
efficients exhibits the same power law mode dependence 
as V: 

Proposition 4 - All the Fourier coefficients f^^f^ of the 
functions fs of the series l^II.12\} satisfy: 



I/, 



is) 



< 



nmkl (^2+^2)3/2' 



for (n, m) ^ (0, 0) . Consequently, for a sufficiently small, 
the Fourier coefficients of the control term f given by 
Eq. nn^) satisfy: 



< 



lnmk\ / 2 1 2^3/9 ' 

for {n,m) ^ (0,0) and for some constant Coo > 0. 

In order to measure the relative magnitude between 



Hamiltonian (|III.16p and /2 or /a, we have numerically 
computed their mean squared values: 



0.13a, 



0.07a^ 



where (/) = dt dx dyf{x, y, t). 



We define the contribution of /2 and /a to the energy 
density by 



62 
63 



1 

8^ 
1 

8^ 



(iv/2n, 
(iv/3n. 



(III.24) 
(III.25) 



For N = 25, these contributions satisfy: 



£2 

£ 



0.1 X 



63 

£ 



0.3 X a^. 



It means that the control terms /2 and /a can be 
considered as small perturbative terms with respect to 
V when a < 1. We notice that even if /a has a smaller 
amplitude than /2, its associated average energy den- 
sity is larger for a of order 1 (more precisely for a > 0.58). 

Remark on the number of modes in V: In Sec lIVI all 
the computations have been performed for a fixed num- 
ber of modes (A^ = 25) in the potential V given by 
Eq. pi.1611 . The question we address in this remark is 
how the results are modified as we increase A^. First we 
notice that the potential and its electric energy density 
are bounded with A' since 



\Vix,y,t)\<a J2 



1 



-,2 oo 



£ < 



8tt 



E 



^ (n^ -|- m^)'^/^ 



< oo. 



< oo. 



(n^ -|- m^)^ 



Concerning the partial control term /2, we see that 
it is in general unbounded with A. From its explicit 
form, one can see that it grows like A^logA^ (see the 
Appendix). Less is known on the control term since 
it is given by a series whose terms are defined by 
recursion. However, from the proof of Proposition O 
in the Appendix, we see that the value a of existence 
of the control term decreases like l/(2A^log A^). This 
divergence of the control term comes from the fact that 
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the Fourier coefficients of the potential V are weakly 
decreasing with the amplitude of the wavenumber. 
Therefore, the exact control term might not exist if we 
increase N keeping a constant. However we will see in 
Sec. IV BI that for practical purposes the Fourier series of 
the control term can be truncated to its first terms (the 
Fourier modes with highest amplitudes). Furthermore 
in the example we considere as well as for any realistic 
situation the value of N is bounded by the resolution 
of the potential. In the case of electrostatic turbulence 
in plasmas kpi ~ 1 determines an upper bound for k, 
where k is the transverse wave vector related to the 
indices n,m and pi the ion Larmor radius,. The physics 
corresponds to the averaging effect introduced by the 
Larmor rotation. 

Remark on the control in the Bohm regime (parame- 
ter a larger than 1): Let us define Vq = V{t = 0) and 
5V = V — Vq. Since the Bohm regime is defined as a 
regime of relatively slow evolution of the potential (with 
characteristic time t^^) compared to the motion of the 
particles (with characteristic time Td), < t^, one can 
introduce as small parameter e ~ t^/ti^ and the inte- 
grable Hamiltonian Hq = Hq-\-Vq, so that H = Ho + eV 
with V = {V — Vb)/e. This approach is only vahd for a 
finite time of order r^^ after which one must redefine Vq, 
Ho and V. 



IV. NUMERICAL INVESTIGATION OF THE 
CONTROL TERM 

With the aid of numerical simulations (see Ref. ^11,] for 
more details on the numerics) , we check the effectiveness 
of the control theory developed in Sec. IIIII by compar- 
ing the dynamics of particles obtained from Hamiltonian 
(|IH.16ll and from the same Hamiltonian with the control 
term /2 given by Eq. (|111.21|l . and with a more refined 
control term /2 + /s where /a is given by Eq. (|III.22|I . 
We use three types of indicators of the dynamics: dif- 
fusion coefficient, Lyapunov indicators, and probabilty 
distribution function (PDF) of step sizes. 

A. Diffusion of test particles 

The effect of the control terms can first be seen from 
a few randomly chosen trajectories. We have plotted 
Poincare sections (which are stroboscopic plots with pe- 
riod 1 of the trajectories of V). On Figures 01 and [3 is 
plotted the Poincare surfaces of section of two trajectories 
issued from generic initial conditions computed without 
and with the control term /2 respectively. Similar pic- 
tures are obtained for many other randomly chosen ini- 
tial conditions. The stabilizing effect of the control term 
(|HL21|I is illustrated by such trajectories. The motion 
remains diffusive but the extension of the phase space 
explored by the trajectory is reduced. 



FIG. 4: Poincare surface of section of a trajectory obtained 
for Hamiltonian 1111.161 using a generic initial condition as- 
suming a = 0.8. 



FIG. 5: Poincare surface of section of a trajectory obtained 
using a generic initial condition as in Fig.l and adding the 
control term 1111. 2H to Hamiltonian 1111.161 with a — 0.8. 



In order to study the diffusion properties of the sys- 
tem, we have considered a set of Ai particles (of order 
1000) uniformly distributed at random in the domain 
< x,y < 1 a.t t = 0. We have computed the mean 
square displacement {r'^{t)) as a function of time 

1 ^ 

i—l 

where Xi(t) — {xi{t),yi{t)) is the position of the z-th 
particle at time t obtained by integrating Eq. ljIII.18|) 
with initial condition X; (0) . On Figure is presented 
(r^(i)) for three different values of a: a = 0.7, a = 0.8 
and a — 0.9. For the range of parameters we consider the 
behavior of {r'^{t)) is always found to be linear in time for 
t large enough. The corresponding diffusion coefficient is 
defined as 

..^lim m. 

t^oo t 

The values of as a function of a with and without 
control term are presented on Figure A significant 
decrease of the diffusion coefficient when the control term 
/2 is added can be readily be observed. As expected, the 
action of the control term gets weaker as a is increased 
towards the strongly chaotic region. 




1000 2000 3000 4000 5000 



t 

FIG. 6: Mean square displacement {r^{t)) versus time t ob- 
tained for Hamiltonian 1111.1611 with three different values of 
a = 0.7,0.8,0.9. 
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FIG. 7: Diffusion coefficient D versus a obtained for Hamil- 
tonian 1IIII.I6I 1 (open squares) and Hamiltonian 1111.1611 plus 
control term l|III.21|l (full circles). 

B. Lyapunov indicator method 

In order to get insight into the action of the control 
term to the dynamics, we apply the Lyapunov Indicator 
method. This method provides local information in 
phase space. It has been introduced to detect ordered 
and chaotic trajectories in the set of initial conditions. It 
associates a finite-time Lyapunov exponent v with an ini- 
tial condition Xq. By looking at the map Xq ^ '^(xq), one 
distinguishes the set of initial conditions leading to regu- 
lar motion associated with a small finite-time Lyapunov 
exponent. The pictures of this map show the phase 
space structures where the motion is trapped and does 
not diffuse throughout phase space, e.g., they highlight 
islands of stability located around elliptic periodic orbits. 

Consider an autonomous flow x — /(x). The Lya- 
punov indicator method is based on the analysis of the 
tangent flow 

^ = Z?/(x)y, (IV.27) 

where Df is the matrix of variations of the flow. The 
Lyapunov indicator is defined as the value 

z.(xo,T)=log||y(T)||, 

at some finite-time T starting with some initial condition 
Xq and a generic vector yo. This definition is very close 
to the one of a finite-time Lyapunov exponent. The plot 
of V versus xq gives a map of the dynamics by highlight- 
ing regions of stability and regions of chaotic dynamics. 
For a chaotic trajectory, the value of the Lyapunov indi- 
cator increases linearly with time, whereas for a regular 
trajectory (periodic or quasi-periodic), it increases like 
\ogt (see rigorous results for nearly perturbed Hamilto- 
nian systems in Ref. [l^)- So in regular regions, this 
Lyapunov indicator is expected to be much lower than in 
chaotic regions. 

Here the Hamiltonian flow is not autonomous. However, 



by considering that t is a new coordinate of the motion 
(and E is its conjugate momentum), we obtain an au- 
tonomous flow with two degrees of freedom. We notice 
that the equations of motion for Hamiltonian l|III.16|) can 
be written as 

X = Re [F(x)e-2^'^*] , (IV.28) 

where 

F(x) = y^V{^, 0) -f- iV-Ly(x, 1/4), 

where ~ {—d/dy , d/dx) and V{'x.,t) is given by 
Ea. HIII.16|) . The non autonomous flow HIV.28|I can be 
mapped into an autonomous flow by considering a third 
equation f = 1. The computation of the tangent flow (of 
dimension 3) follows from the matrix of variations of the 
autonomous flow. We have chosen the third component 
of the vector y following the evolution of the tangent flow 
(|1V.27|I equals to one, which can be done without loss of 
generality since Eq. ljIV.27|) is linear in y. Therefore, it 
reduces to the evolution of a two dimensional vector y 
which is given by 

y = Re [G(x)e-2^^*] y -I- 27rlm [¥{-x)e-'^"'] , 

where G is the two-dimensional matrix G'(x) = _DF 
(matrix of the variations of the vector field F). Fig- 
ure |S1 shows the value of i/(xo, T) as a function of T for 
T G [0, 140] for three initial conditions xq : one strongly 
chaotic Xq ~ (0.865,0.39), one weakly chaotic xq = 
(0.8766,0.39), and one quasi-periodic xq — (0.895,0.39). 
The plot of the Poincare sections of these three trajecto- 
ries up to T = 1000 are shown on Fig. |21 These figures 
show two chaotic trajectories for which there is an over- 
all linear increase of the Lyapunov indicator, and one 
quasi-periodic motion (trapped around an elliptic peri- 
odic orbit). We notice that not only the method is able 
to discriminate early between regular and chaotic mo- 
tions but it is also able to detect weakly versus strongly 
chaotic trajectories for rather small values of T. 

With the control term /2 given by Eq. (|111.21() . the 
equations of motion can be rewritten as 

X = Re [Fe-'*"*] -f f2(x), 

where f2(x) — The equation of evolution of y be- 

comes 

y = Re [G(x)e-2-*] y + 2^Im [F(x)e-2"*] + Dh{^)y. 

With the control term /2(x) -I- f3{x,t) where /a is given 
by Eq. HIII. 221) . the equations of motion can be written 
as 

i = Re [(F + £3)6-2^^*] +f2(x), 

where f3(x) = V-^/3(x, 0) -|-iV-^/3(x, 1/4). The equation 
of evolution of y becomes 

y= Re[(G + g3)e-'*"*]y 

+ 27rlm [(F + {3) g-^-*] + Df^y, 



11 




FIG. 11: Lyapunov indicator map !/(xo,r — 200) for Hamil- 
tonian 1111.1611 for a = 0.4 with the control term /2 given by 

Eq. mrm . 



FIG. 8: Values of the Lyapunov indicator z^(xo, T) as a func- 
tion of time T for three trajectories obtained for Hamilto- 



nian ljIII.16|l with a = 0.4 for the following initial condi- 
tions: one strongly chaotic (a) xo = (0.865,0.39), one weakly 
chaotic (&) Xo = (0.8766,0.39), and one quasi-periodic (c) 
Xo = (0.895,0.39). 
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FIG. 9: Poincare sections of the three trajectories of Fig. |H| 



where 33 (x) = Dij. 

On a grid of 10000 initial conditions Xq S [0,1]^, we 
compute z^(xo,T) for T = 200. Figure [TUl represents 
the Lyapunov indicator map for a = 0.4 without control 
term. Figure ITTI shows for the same values of parameters, 
the Lyapunov indicator map with the control term /2. 
The general effect of the control term /2 is a decrease 
of the magnitude of the Lyapunov indicators. However, 
the stabilization effect of the control term /2 is not uni- 
form. There are regions where the (partial) control term 
/2 fails to stabilize the trajectories, e.g., in the region 
near Xg = (0.1,0.4). Figure [T21 plots the Poincare sec- 
tion of the trajectories starting at Xq — (0.085,0.385) for 
Hamiltonian (|III.16|I with and without the control term 
/2 given by Eq. I|III.21(I . Wc notice that this trajectory 
is more chaotic and more diffusive with the control term 
than without. 




FIG. 12: Poincare sections of the trajectories of (a) Hamil- 
tonian I ini.l6l l and (b) Hamiltonian 1111. 16II with control 
term /2 given by Eq. 1111. 21L with initial conditions xo — 
(0.085,0.385). 



In order to see the global stabilization effect of the par- 
tial control term /2 , we notice from the values plotted in 
Fig. ^1 that for a = 0.4 about 25% of the trajectories 
have a Lyapunov indicator less than 5 at T = 200 with- 
out the control term /2, compared with 70% with the 
control term (from Fig. 111(1 . Figure El represents the 
histograms of the Lyapunov indicator at T = 200 for 
a = 0.4 with and without control term. The first peak 
in the upper and lower panels corresponds to the regular 
component of the phase space. We clearly see that the 
second peak corresponding to the chaotic component is 
drastically reduced with the addition of the partial con- 
trol term. The effect of a more refined control term is 
observed in Fig. 1141 for a = 0.6 when the control term /2 
starts to fail to reduce significantly the chaotic part of 
phase space. We see that the proportion of the regular 
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FIG. 10: Lyapunov indicator map v{-x.q,T — 200) for Hamil- 
tonian IIII.I6II for a = 0.4. 



FIG. 13: Histograms of the Lyapunov indicators at T = 200 
for a — 0.4 with and without control term computed in 
Figs. ITHl and im 
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FIG. 14: Histograms of the Lyapunov indicators at T = 200 
for a — 0.6 without control term (upper panel), with control 
term /2 (middle panel) and with control term + /a (lower 
panel) . 



FIG. 15: PDF of the magnitude of the horizontal step size 
for Hamiltonian 1111.161 with a — 0.7 without the control 
term (open squares), with the control term /2 (fuU circles), 
and with a truncated control term with twelve modes (open 
circles) . 



trajectories has increased with the addition of f^. More 
quantitatively, 8% of the trajectories of the Hamiltonian 
without control have a Lyapunov indicator smaller than 
7 at r = 200. With the addition of /2, this proportion 
is increased to 25% whereas it is around 30% with the 
addition of /a. 



C. Horizontal step sizes 

In order to investigate the effect of the control term on 
the transport properties and its relationship with single 
trajectories, we have computed the Probability Distribu- 
tion Function (PDF) of the step sizes. Let us define the 
horizontal step size (resp. vertical step size) as the dis- 
tance covered by the test particle between two successive 
sign reversals of the horizontal (resp. vertical) component 
of the drift velocity. The effect of the control is analyzed 
in terms of the PDF of step sizes. Following test particle 
trajectories for a large number of initial conditions, with 
and without control, leads to the PDFs plotted in Fig. 1 151 
for Hamiltonian l|III.16|l without and with control term 
HIIL2ip for a = 0.7. A marked reduction of the PDF is 
observed at large step sizes with control relatively to the 
uncontrolled case. Conversely, an increase is found for 
the smaller step sizes. The control quenches the large 
steps (typically larger than 0.5 for a = 0.7). Such a re- 
duction of the probability to achieve a large radial step 
will modify the transport efhciency and in particular re- 
duce the diffusion coefficient. 

In order to give support that the first peak of the his- 
togram of the Lyapunov indicator (see Fig. I13|l is asso- 
ciated with the small step sizes, we have plotted on the 
Fig-Elthe distribution of horizontal step sizes of the tra- 
jectories with a small Lyapunov indicator (smaller than 
7), and also the same PDF for trajectories associated 
with large Lyapunov indicator (larger than 7). This re- 
sult gives support to the fact that the control term re- 
duces the diffusion of trajectories by reducing chaos in 




step size 

FIG. 16: PDF of the horizontal step sizes of the trajec- 
tories with small Lyapunov indicator (smaller than 7, full 
squares) and with large Lyapunov indicator (larger than 7, 
open squares) for Hamiltonian Hill. 1611 with a = 0.4. 

the system (by the creation of invariant tori [2^ , see also 
the global picture in the conclusion). 

V. ROBUSTNESS OF THE CONTROL 

In the previous section, we have seen that a truncation 
of the series defining the control term by considering the 
first or the two first terms in the perturbation series in 
e, gives a very efficient control on the chaotic dynamics 
of the system. 

In this section we show that it is possible to use an 
approximate control or to make a small error while 
computing the control term and still get an efficient 
control of the dynamics. 

Below, we give numerical evidence for the following 
statements : the reduction of the amplitude shows that 
one can inject less energy to achieve a significant control. 
The truncation of the Fourier series indicate that one can 
simplify the control term and still get a significant con- 
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FIG. 17: Diflusion coefficient D versus the magnitude of the 
control term 1111.2111 for a = 0.7. The horizontal dashed line 
corresponds to the value of D without control term. The 
dash-dotted line is a piecewise linear interpolation. 



trol. The change of phases shows that one can introduce 
some error in the phases and still get a significant control. 



A. Reduction of the amplitude of the control term 

We check the robustness of the control by increasing 
or reducing the amplitude of the control 24] . We replace 
/2 by (5 • /2 and we vary the parameter 5 away from its 
reference value 5—1. Figure [T7I shows that both the in- 
crease and the reduction of the magnitude of the control 
term (which is proportional to 6 ■ o?) result in a loss of 
efficiency in reducing the diffusion coefficient. The fact 
that a larger perturbation term - with respect to the 
computed one - does not work better, also means that 
the control is "smart" and that it is not a "brute force" 
effect. 

The interesting result is that one can significantly reduce 
the amplitude of the control {5 < 1) and still get a reduc- 
tion of the chaotic diffusion. We notice that the average 
energy density 62 ((5) associated with a control term 6 ■ ji 
is equal to ei[S) =5^62, where 62 is given by Eq. (|III.24p . 
Therefore, for 5 = 0.5 where the energy necessary for the 
control is one fourth of the optimal control, the diffusion 
coefficient is significantly smaller than in the uncontrolled 
case (nearly factor 3). 



FIG. 18: Contour plot of the truncation of j^ix, y) for a = 0.7 
containing the twelve Fourier modes of highest amplitude. 



where fnX is 



(2) 
nrn 



{nnii — nim)e 



3/2 



a x - 

^ nh^, K + ml)y^ [{ni - nf -f (mi - mf)] 

(V.30) 

The truncation of the Fourier series is made by consider- 
ing the Fourier modes with an amplitude greater than or 
equal to e, that is the sum in Eq. (|V.30|) is restricted to 
the set of modes (n, m) such that \fnm\ > £• For exam- 
ple, if we consider e ~ 9 • 10"'', there are twelve modes 
in the sum, which are the modes with wave vector (0, 1), 
(0,2), (1,-1), (1,-2), (1,0), (2,0) and the opposite wave 
vectors (/2 is real), compared with the total number of 
modes of the full /2 which is about 2000. Figure ITSl 
shows the contour plot of /2 obtained with only these 
twelve modes. Figure 1191 shows the diffusion coefficient 
for the dynamics of the truncated control term versus 
the number of Fourier modes kept in the truncation of 
/2. The reduction of the diffusion (with respect to the 
uncontrolled case) also holds for a very simplified con- 
trol term containing only the few highest Fourier modes 
of the full control term. If we replace /2 by the trunca- 
tion with twelve modes we see that the effect is still a 
strong reduction of the diffusion coefficient, a reduction 
of about 25%. The energy density 62"^^'' (see Sec. IIII C|) 
of this truncated control term with respect to the energy 
density £ is 



0.009 X 



(V.31) 



B. Truncation of the Fourier series of the control 
term 

We show that a reduction of the number of Fourier 
modes of /2 can still significantly reduce chaotic diffusion. 
The Fourier expansion of the control term /2 given by 
Ea. ()iii.^l|) is 



/2=E/' 



(2) 2i7T(nx-\-my) 



(V.29) 



where £ is given by Eq. (|111.23|l . that is less than 1% of 
the energy associated to the electric potential. It is inter- 
esting to notice that the energy density of this truncation 
with respect to the one of the full control term /2 is 



J12) 



62 



0.09, 



(V.32) 



where €2 is given by Eq. Ijlll.24|l . that is less than 10% 
of the energy associated with the full control term /2. 
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FIG. 19: DifFusion coefficient D versus the number of Fourier 
modes n in the truncation of the control term /2 for a = 0.7. 
The dashed line corresponds to the case without control term, 
the solid line corresponds to the value of the diffusion with 
the full control term /2 and the dash-dotted line corresponds 



to a power law interpolation (oc ' 



-1/2 



FIG. 20: Diffusion coefficient D versus the phase error 7 in 
the expression 1111. 2H of /2 for a = 0.7. The dashed line 
corresponds to the case without control term, the solid line 
corresponds to the case with the full control term and the 
dash-dotted line corresponds to a quadratic interpolation. 



Moreover, we see from Fig. 1151 where the PDF of hori- 
zontal step sizes is plotted, that the effect of this coarse 
grained control term /2 reduced to twelve modes is also 
to quench large step sizes. 

More generally, these results show that a partial knowl- 
edge of the potential, e.g. on a grid (coarse grained) is 
sufficient to obtain a significant control of the dynamics. 

C. Change of phases in /2 

We check the robustness of the control with respect 
to an error introduced in the phases of the control term 
given by Eq. Ijlll.2ip . i.e. we change the phases by 
(^„™ in Eq. (|III.2ip by : 

<^«m = ^7iva + 7 • "^nm > (V.33) 

where f'^^ are uniformly random distributed phases in 
[0, 27r] , 7 is the amplitude of the error and ^nm are the 
correct phases. Figure [501 shows the diffusion coefficient 
versus the phase error 7 for a fixed value of a (a = 0.7). 
We notice that the chaotic diffusion is still significantly 
reduced by the control with a small error on the phases. 
The diffusion coefficient is still strongly reduced by a fac- 
tor greater than 2 for a phase error of 5%. For small 
values of 7 the diffusion coefficient versus 7 is well fitted 
by a quadratic interpolation, that is I?(7) = Dq + Di^^ . 

VI. CONCLUSIONS 

We have provided an effective new strategy to con- 
trol the chaotic diffusion in Hamiltonian dynamics using 
small perturbations. Since the formula of the control 
term is explicit, we are able to compare the dynamics 
without and with control. The idea of the control is pic- 
torially represented in Fig. [2] A Hamiltonian Hq + eV is 



controlled by adding a control term /. The naive choice 
for a control term would be / = —eV but this would be 
useless since it is of the same magnitude of the source of 
chaotic transport and thus would require a major modifi- 
cation of the physical condition of the system of interest. 
In this article, we have presented a way to design an in- 
tegrable controlled Hamiltonian He with a small control 
term / of order e^. This controlled Hamiltonian is con- 
jugate to Hq (we assume for simplicity that TZV ~ 0). 
This construction of the controlled Hamiltonian works 
well up to some value ei. Moreover, we have shown that 
the control is robust, in the sense that one can use an 
approximate controlled Hamiltonian H,. which is not in- 
tegrable but -being sufficiently close to He - generates a 
more regular dynamics (presence of invariant tori) with 
respect to i?o + ■ For instance, we have shown that 
one can successfully use a truncated control term of or- 
der and that one is allowed to tailor it to some specific 
requirements on its shape, on the energy necessary to 
achieve control and also according to a partial knowl- 
edge of V{e.g. a truncation of the Fourier series and an 
error on the phases). 

The invariant tori that have been created by adding the 
control term / of order are those which were broken 
by increasing the amplitude of the perturbation, mean- 
ing that these tori are those of a Hamiltonian Hq + e'V 
where e' < s (up to some smooth canonical transform 
close to the identity transform). In order to illustrate 
this statement, we have plotted in gray two regions of 
existence of a given invariant torus (specified e.g. by 
its frequency) . The uncontrolled Hamiltonian Hq + eV 
does not have this invariant torus whereas the controlled 
one He — Hq + eV + e^f does. The controlled Hamil- 
tonian He is conjugate to the controlled Hamiltonian 
H'e = Hq + s'V + e'^f for e' < e (since they are both con- 
jugate to Hq). Since Hq + e'V is inside the ball around 
the integrable Hamiltonian H'e, the invariant torus of 
Hq+e'V of the selected frequency is a small deformation 
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FIG. 21: Global picture of the control: The bold curved 
segment of curve represents a set of integrable Hamiltonians 
around Ho- The gray circles are the domains of existence of a 
given invariant torus around an integrable Hamiltonian. The 
gray arrows represent two ways of controlling the Hamilto- 
nian Ho + eV: the first one of order e and the second one of 
order e^. The arrows from He to H'^, and from H'c to Ho + e'V 
represent close to identity canonical transformations. 



of the torus of the controlled Hamiltonian H'^ and hence 
a small deformation of the torus of the controlled Hamil- 
tonian He- Therefore the invariant tori of He obtained 
by means the control are small deformations of the tori 
of Hq + eV which were broken by increasing e. 

We have applied this general technique of control to a 
specific model, describing anomalous electric transport in 
magnetized plasmas. In particular, we have shown that 
the control term is robust, meaning that one is able to 
simplify it, to reduce its amplitude or to make a small er- 
ror without changing its overall action of reducing chaotic 
transport. Even though we use a rather simplified model 
to describe chaotic transport of charged particles in fu- 
sion plasmas, our result makes us believe that through 
some small smart modification of the electric potential a 
relevant reduction of the turbulent losses of energy and 
particles in tokamaks could be attained, for the moment 
at least in principle. 
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APPENDIX A: PROOF OF PROPOSITION |31 
ANDH 

1. Proof of proposition |S| 



The terms of the series l|11.12|) can be written in the 
following form: 



is) 



27r(nx -I- my) + (Pj^l^^ - 27rfci 



where the sum over k is from — s to s, and the two sums 
over n and m are from —sN to sN. From the recursion 
formula (pT3)l . we have 



N 



mn — nm 



Jn'ni'k' 2s ^ (n2+m2)3/2 

n.m— 1 ^ ' 

I •> — n,m' —m^k' — 1 n— n',m— m',A;' + l j ' 

(A.l) 



We use the following norm : 



\fs\\ - sup \fi'l,\ 

n,m,k 



From Eq. ljA.l|) . we get 



N 



\fs\\ < sup > 

S I t ' ' 



\mn' — nm'\ 

n',m' ^ , (n^ + to2)3/2 
' n,m— 1 ^ ' 

Since, \mn' — nm'\ < sN{m + n), we have 
ll/.ll <aA||/._i||, 

where 

JV 



A = 2iv y — — — 

n,m— 1 ^ ' 



(A.2) 



It follows that 

IIMI < (aA)^-ia2-3/2, 

since ||/i|| = a2-3/2. Therefore, the series (|1L12|I con- 
verges for a < 1/A. We notice that for N — 25, A « 135, 
i.e. the series converges for a < 7 x lO^'^. 
Since n i— > m/ ['n? + m?)^/"^ is a positive and decreasing 
function, we have 



mdn 



< 



mdn 



N AT 

Em I 
[V? + m2)3/2 - y„ (^2 to2)3/2 (^2 + ^2)3/2 

By rescaling the integral (t = n/m) and using the fact 
that ^1)3/2 = 1, we have 



A < 27V V — < 27VlogiV-|-27iV+ 1, 



(A.3) 



where 7 is the Euler-Mascheroni constant. In particu- 
lar, we notice that the bound on A increases like iVlogiV. 
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2. Proof of proposition 

Concerning the regularity of the functions fs and /, 
we would like to show that 



(A.4) 



We notice that Eq. 1A.4|I is satisfied for fi^V given by 
Eq. HIII.16|) for C > 1. The inequahty gives 



\f^l,'\<-'- 



N 

E 



|mn' — nm'\ 



^ {n? + m2)3/2[(n' - n)2 + (m' - m)2]3/2 ' 



We can always assume that both n' and m' are positive 
since we have |n' — n| > ||n'| — n| which gives 



N 



m + n 



^ (ri2 + ^2-)3/2J(-|j^/| _ „)2 _|_ (l^^/l _ ,7j)2]3/2 ' 

where we have used the inequality \'mn' — nm'\ < sN{'m+ 
n). 

We have to distinguish the following cases: (i) n' > N 
and m' > N, (ii) n' < N and m' > N, {in) n' < N and 
m' < N . We notice that by symmetry the case n' > N 
and m' < iV is similar to (ii). 

For n' , m' > N, we have n' < {N + l)(n' — n) and m' < 
{N + l){m' — m). Thus we have 



< 



[n' — n)2 + (to' — to)2 ~ n'2 + to'^ ' 
which leads to 

(,) I ^ a^C^-i(jV+i)3A 



I n'm'k' 



(n'2+m'2)3/2 ' 



where A is given by Eq. ljA.2|) . 

For n' < N and to' > TV, we have to' < (A^ + l)(m' — to). 
By using the estimate (n' — n)^ + (m' — to)^ > (to' — to)^, 
we have 



1 



< 



(iV + l)= 



(n' — n)2 + (to' — m)2 to 
and since n' < to', 



(A^+ 1)2 ^ 2(A^+ 1)2 



/2 



71^ + TO ■ 



Thus we have 

(,) a-C-i23/2(Ar + l)3A 

l-^"'m'fe'l- (^/2+^/2)3/2 

For n' < N and m' < N, we use the crude estimates 



(n' — n)2 + (to' — m)^ 



< 1, 



and 



Thus we have 



1 < 



2N^ 



n'2 + to'2 



(,) ^ a-C-i23/2jV3A 

l.'n'm'fc'l - (^,2 _^^/2)3/2 " 



By denoting C = 2^^'^{N + 1)^X, Eq. ifOj) is satisfied for 
all n, TO and for all s. 

It follows that for a < 1/C, the same inequality holds for 



(n2 + to2)3/2 



where Coo - a2CV(l - aC). 
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